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CALOULATIOH OP POTEUTIAL PLOW PAST AIRSHIP BODIES 111 YAW* 

By I. Lot z 



The calculation of tlie potential flow past an airship 
in yaw is accomplished, according to Von Zarman, by divid- 
ing the air stream into a flow parallel to the airship ax- 
is and one perpendicular to it, after which the potential 
flow about the solid of rotation exposed in axial direc- 
tion can then be determined by a method suggested by Ran- 
kine, The"~^is of the body is superposed by sources and 
sinks and the intensity so defined that the closed stream 
surface formed by superposition of parallel flow on this 
source- sink flow is coincident with the surface of the 
solid. Von Karman (reference 1) pursued an analogous meth- 
od for defining the potential set up by the lateral flow. 
He disposes dipolars** on the axis, so that the intensity 
of this superposition is again defined by the postulate 
that for the fluid^otion, produced by superposition of 
parallel and dfpoTar flow, the surface of the solid is 
stream surface* Both problems yield int.egral equations. 
For^^^}proximate solution, he substituted a definite number 
of pie'^es with constant superposition for the continuous 
superposition. A numerical calculation of examples accord- 
ing to this method revealed certain difficulties for sol- 
ids in yaw with not very slender bows, which finally prompt- 
ed tho use of sources and sinks on the surface of the solid 
instead of dipolars on the axis, and to compute the field 
of flow by this a j f&ntfy . Tho experiences collected from tho 
calculations and which should prove useful for othor sim- 
ilar problems, form the subject of this report. 



tt 

*Zur Berechnung dor Potential stromung um quargestellte 
Luf t schlff korper-, .Ingenieurr^Archi-vf , Vol, II, 1931, pp. 
507-527. 

**0n Dipolars and Their Fields of Flow, See Hutte, Vol, I, 
26th edition, p. 364, 



N.A.C.A, Technical Memoranduni lTp« 675 

I. AXIAL SUPERPOSITION 
A. Partially Constant Dipolar Superposition 



The axis of the airship is assumed to coincide witli 
axis X, counted from the ■bow. Cylindrical coordinates 
are introduced so that 

z = r coBcp , y = r sin<p (Pig. 1.) 

The solid itself, of length I is assumed given hy its me- 
ridian .curve r(x). 

The normal velocities induced on the surface oy the 
parallel flow inairection of^liegative axis z, are pro- 
portionate to cos cp. Since dipols lying on the ari s of 
the solid also set v.p on the surface of the solid normal 
velocities proportionate to cos q),-l*srangle 9 does not ap- 
pear in the integral equation defining dipolar motion \i. 
This equation rests, as already stated, on the premise that 
the surface Is^stream surf ace, so that on it the normal 
velocity, composed of^ipolar motion and'^tiomp onent set 
hy the flow, disappears. Designating the point of the ar- 

is carrying the superposition "by 5*, expressing ~ ' 

and W = velocity of air flow, this equation for the ad- 
herence of flow to a parallel of latitude with abscissa x 
reads as 

. / .... i ^(e) M . ^ + 3 r- / tiCllA^ , 

° /^'+ (x - ° + (x - i)^" 

, - 3 r r» / — _ -=rrr r = 4 tt W (l) 

o -y/rs + (x - g) 

It is known, and Earman also pointed it out in his 
summary, that exact substitution of the "body "by a superpo- 
sition on. the axis of rotation is possible only "when the 
analytical continuation of the potential function free of 
singularities in the outer space of the solid can he ex- 
tended to the axis of symmetry without encountering singu- 



*The zero point of | coincides with that for x. 
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^v^SiT distribution the 



lar points," Since the ^''J-estiori of when or ^rhen not this 
continuance is pos.si"ble, has nisver Iseen answered, nxtmer- 
ical calculations are apt to have di sagreeatle surprises, 
in store:.*. ■• 

Fith n -epste*^ of constant 
flow adhe.res. only in n points (control points) of the 
meridian curve; that means on n circles of the hody pe- 
riphery. (Fig. 2.) The Tesult is an n system of liuaar 
■equations for the n ■, unknown dipolar intensities tax?.. 
With the: vth dipol space reaching from. Iu to tv+i 
(fig. 3) , the control point x (reference l) for flow** 
adherence is expressed as . . 



2u M-u [ A Su + r« A fu] = 4 tt r^ W 



::(2) 



wherehy 



where 



Agy = ei^*iv) - sCx.Su+O*! 

A f v = f(x,lu) - f(3C,lu+i),J 
g(x, I) = cos^ (2 - cos2 ^J-),^ 
f (x., I) = sin3 ^ 



(3) 



(4) 



wi th 



.tan ^ = 



X - i 



(yig.;. 3.) 



Karman has computed the functions g and f and' the; 

X - I 



appear as dependent on 
equation (2) "becomes 



For point x =.9, r = Q 



2 1 Ml 



v+i 

where m = limit (r r*) 
X— >0 
r S 



m 



P W (-~-- -^)=. W- (5) 



^1? 



'♦In the •Jablonowski Prize Essay, 1914, G-, HerglOtz treaited 
a similar pro'blemi Analytical continuanee -of potential into 
the inside of the adhering masses, although the fundamental 
difference of the pro'blems makes it impossi'fale. to expect 
any explanation, 

♦♦Manifested in (l) by inserting partially constant super- 
position 



4 
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m is a finite number, /because, r ordinarily goes, to zero 
as c^/~x'. Then the product* (r r«) for small x is 



that is, finite for x — »-0 and equal to the radius-, of the 
curved sphere for point .x = 0, r = 0, 

A different method of visualizing the dipolar inten- 
sities is to cut the solid parallel to its axis and perpen- 
dicular to the parallel flow and stipulate that the flow 
through the outside of this section, generated by the di- 
polar motion he equal to the flow of the undisturbed stream 
through this area, i«e«, equivalent to the. section multi- 
plied by the flow velocity at infinity. It is not advisa- 
ble to effect n such cuts to define the n dipolar in- 
tensities, because the slender forms of interest to us 
render the exact graphic analysis of the section difficult 
and at the same time becloud the dipolar effect when the 
-cut is made at a di,stance r contiguous to Tmai ^rom 
- the axis. But it is advantageous to. occasionally replace 
one or two equations with such flow conditions in the sys- 
tem (2). 

The length of the constant dipolar intensities must 
be decided at the very beginning. Knowing |j, to be about 
proporti^^al^o rs according to Von^^arman, the length 
of the s paces •-depenas on the va - 4M>'n-s-i-^by of -^te change of r. 
(j^'' will be less toward bow and stern. The major obstacle 
lies in the choice of the starting point of the superposi- 
tion, whose distance from the bow may be. assumed x^ , For 
the sphere the ratio of Xj^ to curvature radius m is = 
1» Long extended bodies having on the axis a stronger di- 
polar motion, the interference velocities outside of the 
body increase conformably to the flow condition, which 
means a shortening of distance x^^ , Consequently, the 
starting po,int must be ahead of the center of curvature 
of the bow; for the rest, its selection must be left to 
experiment. 

If a calculation is made first it a.ffords a check on 
whether or not point x^^ has been correctly chosen. If 
the surface of the body is reallytJ surf ace of flow then 
even the higher differential quotients r", etc., of the 



the' derivat-i<wi of the 



(r r»)=-^-^/; that is, . r r» = 
dx 

curve TT r2 (x) divided by 2 tt. 
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meridian' section of the bio sed ■ stream surf ace-j prodtLced 
Tjy superposed dipolar and parallel flow, must correspond 
with those of the surface of the ijody, ITe shall see wheth- 
er this igt>ie. case "for the bow. Differentiate equation 
(l) MooPv^tni^ t6 3c . and form' the limit for x = Q and 
r = Q . Putting ■ • • ■ ■ • • 



and 



limit [ (r'«'2 + r r") 1] = 3' 
x-*0 



r — >0 

it manifests 

p I^t; |se3u + 7 m.e^y- 5 65^ != 0 (9) 

The condition is extremely sharp. It shows that even 
a minor error in the selection of is immediately fol- 

lowed Toy a tangible deviation from zero at the right side. 
Even an estimate of the anticipated ]x values may, under 
certain circumstances afford a check on the ^fulfillment of 
this condition,. -The Gp-^, value s drop very .' quickly , so 
that for the balance onl.y the first spaces come into, con- 
sideration. 

The speedy abatement of tho 9pu valuos -with rising 
p induces the :;.f ollowing : ■ Assuming that very high •d**#sn» is4vw«^-,.. 
IfcioavB.. of -the bulkhead surface curve are known, it '.would'. be 
possible to compute directfi-jby repeat*»§^ differentia- 

tion- .of equation (l) and formation of its limit, for x .= 0, 
r = 0, after having proceeded far enough so that . Q^y. 
abates rapidly, so that 11^, the first term of the sum,, only, 
need be considered, 

Uow the course of the dipolar motion for pointed 
stern must bo included yet. How close to the tip i3 the 
superposition to extend? It is readily shown that inte- 
gral equation (l) 



TO 



,y r8 + (x - e)» 0 y r» + (x - o 



. 3 r r« / ^x(e) (X - l)__d i_ ^ 4 „ 
.0 Jl^+U-if 



6 
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which the dipol motion must fulfill, affords no positive 
solution at the right of the maximum section for r r' = 0 
vrhich 'begins in s = I, r = 0 as |J,(i) = c(l - ir^ , 
where X and c are positive figures. It means: the di- 
polar motion cannot reach the tip and dipolar motion re- 
sults only in an approximated rounded-off sterra. In a 
practical sense there is, of course, no exact tip. Taut 
merely a point of very high, hut still finite curvature. 
The question of final point of the superposition then 
would, provided a certain roundness is desired, revert to 
the considerations made for the "bow. But the superposi- 
tion toward the stern will he so small that a minor vari- 
ation in final point will have no effect on the total re- 
sult (on the moment, for instance). 



3, Continuous, Partially linearly Variahle 
Dipolar Superposition 



Jor not very slender solids the dipolar superposi- 
tions defined, according to the descrihed method are very 
much contingent upon the chosen space division. (See figs, 
5-8.) This fluctuation admittedly diminishes .as the num- 
"ber of dipolar spaces increases. But for mathematical 
reasons a minimum number is desirable, in ■ order to keep 
the order of the resolvable equation system low. These 
considerations were the reasons for assaming a continuo^is, 
partially linearly variable superposition in place of a 
partially constant dipolar motion. The result was, as 
proved by examples to be given elsewhere, quite satisfying 

Proceeding from the general equation (l) for flow ad- 
herence, we first develop the equation system. With 

p. (I) = a-v + ^ (10) 

in the I'th interval, which extends from to ^-p-l-i 

(fig, 4) , equation (,l) becomes 

t 



- 2V / ^ d|+2r^ Zi? / ~™~:z:-5- ^ ^ 

^ iv -f r2 + (x^ ' ^^2 + (x-e)^ 



1 /r2+(x-e)^ 
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or , when writing " 

B-v + >v I t (au + b^ x) - hu. (x - I) 



2U (ap+Ta^jX) 
1 



-3rr» / 



n 



=4tt¥ (11) 



The integrals of the first sum on the left side appear 
by Von Karman on page 3; those of the second sum are read- 
ily evaluated. Multiplying the whole equation by. r^, 
yields 

11 n 

(au+bpx) [A gy+r' Af^D " by r [ - (sin^+i -sin ^y) + 
1 1 

+ (sins^+i-sins^u) +r ' (cos3^u+i-cos3^v)] = 4TTr2¥ (12) 
For abbreviation, we write; 

A 1^ = (sin3^- sin3i9x)+i) -(sin^ - sin^+i) (13) 

A ~ (coss^u- cos^^^u+i) t 

so that (12) becomes 

n a 

T.V (a^+byx) [Agu+i" Afp] + T.^^ b^r [A +r ' A kp] =4'rTrSW . 

(14) 



1 1 
Carrying the simplification further. 



A Sy. + r ' A f^ = , 'I 
x(Agi, + r« Afp) + r (A Iv + r« Ak^) = 11^; J 



(15) 
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Equation (14) , which defines and- "b^, assumes the- 

form 



Zv (au Ip'+'Tju 11^) = 4 TT r2 W 



(16) 



The superposition begins in with- value c. Then 



ai = c - ^1 



= c + bids - ii) - Is. 



> (17) 



au = c+bi (Is- Ii5+b2 (Is- Is)^' • • •+^u.-i ( ly- lu-i) -"bulu J 
From (16) follows 

+^2 . {(ia-U) 2V I^- Is I3+ lis} 



• + 



(18) 



+^n- 1 { ( In- ^n- 1 ) In- ^n- 1 In- 1 +1 ^n. 1} 

+^ { - ^ In +IIn}= 4 TT r2 W. 

"which is the equation system for computing the (n + 1) 
unknown factors ■bj.,,.bj^. For the point x = 0, r = 0, 

the relation with the notations of equation (8) is 



|u |a^,(egu-3me3^)+bp(eip-3me3T;)| =4ttW 



Putting 

6 gp - 3 m Ggu = IILp, 
6.1U - 3 m = IVb . j 

reveals for point x = 0, r = 0 



I 



(19) 
(20) 
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- n . . . n ■ ; ■ 

c 2y 1.1 iu +'b'^ ^ ( ^ -\ ) zu'-xi lu - i i ' i ^'h } "' ' 

{(1.3.-18.) Illy- Is Ills +1^} 
• + ■ 

+tn ■'*■' .•' ; -in ' l^ln +IVn-j. =4TrW (21) 

It should "be observed', when'' carrying out the calcula'- 
tion, that the functions. . g,. f , k, I are the same for all 

bodies dependent on s— ^— ^; or, in other words, construe 

k and I in addition to g and f as in the Karman method, 
The , development of equations (18) and. (21) takes somewhat 
longer than when assuming partially constant, dipolar su- 
perposition, but the decidedly better results- by less num- 
ber of spaces compensates for this drawback. 



-.'jC, Examples for. A and B. 



.To illustrate the preceding expositions there follows 
a compilation of calculated dipolar superposition for bod- 
ies which in axial flow are to be substituted by simple 
singularities on their axis. The body contour is indicat- 
ed and the- sj-jraputed dipolar superposition plotted in . units 
of the supc'-r-position, which corre spends to an indefinitely 
long cylxn-ipc of diameter (2'E;). of maximv.m cross section 
in a transverse flow. Designating the plotted values with 
[i,^ (i) , the dip.ol moment is given by 

,H(|) = 2 TT R» W (E) 

(1(7 == rate- of flow) i The abscissas of the utilized control 
points are shown by X-la. 



*The functions g, f, k, I were omitted here because of. 
the limited space for the scale, and the comparative ease 
with which the functions can be plotted. 



10 
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1, Figure 5.: a "body which is produced lay. a source 
and sink of fj.'bsolutoly equal yield in a flow parallol to 
its connecting lino. The yield in relation to the axial 
rate of flow, is so p.roportipned that the maximuia diameter 
equals' 1/5 of the distance of the singularities. Bow and 
stern have the -s-am^" -sh-ape-, so the how is shown only. It 
is readily apparent from Pigure 5 that the starting point 
of the superposition: was chosen too far away from the "bow. 
As a result "of its too remote distance from the how the 
superposition of the first spacing is very large and, in 
order to neutralize its effects on the other parts of the 
contour., a. negative superposition must ensue, 

2. ITigure 6! hody produced as in Figure 5; yield of " 
source and sink such as to' make the maximum diameter equal 
the spacing of the singularities, 

a) The numher of, dipolar spacings is insufficient, 

b) and c) 'Partig,lly -linearly yariable. .st eady dipolar 
superposition assumed. In case c) it stipulates compli- 
ance with differentiated limiting conditions (compare page 
5), Note that the fluctuations of the superposition in. 
a) are perfectly in accord with rate of ' superposition 

c).' 

, 3, Figure 7:. Fuhrmann (reference 2) shape Ho. IV. 
(The sketch shows the source distribution producing the 
body in axial flow.) The minor negative dipolar superpo- 
sition adjacent to the bow is most probably due to not 
q^tite exact selection .of starting point of superposition, 

\a^*4, Fuhrmann shape lTo« 1, This body has a very blunt 
nose. The dipolar superposition is computed in three 
ways. Because of the slight effect of nose and stern in 
very slender bodies the dipol superposition can, so to 
speak, be effected for two semi-bodies*. For this body 
the dipol superposition, cojaputed according to Von Karman ' s 
method, manifests unusual fluctuations at the bow (a). 
For that reason the equations, given by the postulated ad- 
herence of flow in the first two control points, were re- 
placed by two flow conditions, which -abated the fluctua- 
tions very materially (b) . In Figure 8c the starting 
point is much farther away from the nose which, being very 



*The 



minor 



mutual effect can be accounted for by 



iteration. 
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blunt, causes the same phenomena as in e:?:ainple 2, 

.As- concerns bodies with blunt nose, the iteration of 
the partial equation system for , the nose is impossible. 

II. SUPBHPOSITIOIT OH THE SUHFACE 

A., .Development , of Integral Equation for 
Surface Superposition and Its Reduction to a 
Unldlmensional Integral Equation 



The hesitation voiced at the beginning of this report 
against substituting for the body a dipol superposition on 
the axis, added to the experience gained by computing the 
cited examples, induced me to apply -a different method for 
defining thei potential $ of the transverse flow around a 
solid of rotation. Using the same notations as before, 
namely, x = axis of body, z = negative direction of par-' 
allel'flow, which' is perpendicular to x, ajid y = direc- 
tion perpendicular to x and z , •^iUpotential 0 is defined 
as follows: It must satisfy the potential equation • A'J = 
0, it must act at infinity as (- W, z) , and on the sur- 
face of the' hody must be = 0 . (n = direction of the 
outside normal). Then resolve $ into $i and $s , where 

$1 = - W z. On the surface ^ " * -"^^^ fulfills 

the conditions for $3 must disappear at infinity and 

the normal d orivati on on the surface is — = + W r; — • 

on rjn 

The problem thus resolves to defining $2. This potential 
was approximated in the Rankine- Von Earman method by the 
dipol superposition of the body. axis. But, conformably to 
the theory of mass potentials j $2 can equally be repre- 
sented as potential of an eleiaentary source- sink superpo- 
sition f spread over the surface of tha-bo^^j^ On the 
precise that the body everywhere has a l^^^jS^BT^^ cLe fined 
tangential plane, that- i s , no - cusps , this superposition 
can be, expressed by (ref.erence 5 ) : .lGnMia::l ) 
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Here S, ..surface of tody. , 

dCJ, ■ surface diffdrential at point |, ..ri,. , 

Y, angle' between connecting line of fixed sur- 
face point P(x,y,z) and running point 
II(i*r\,0 with, the normal in point P(x,y,,2), 
(Compare.'f ig, 9,) 

R, distance of fixed from running point. 

The "body "beiing symmetrical in rotation, cylindrical 
coordinates are introduced; 



fixed point P: r= y z2 +y3 , z=r cos Cfjp, x, 

running- point TT: P= J~^^?rf , £-P cos q)p, i, (23) 

so that (22) becomes,. 

^ cos i[f W cos <p« 

2iTf (x,r,q)j.)+ / f (|,P,9p) -rri <i<> = ,. = (24) 

■ , /^W • 

dr dp 

In addition it is written - ~ Extrapolation 

then yields from the geometrical relations shown in Figure 
9, . 

[r-r'(x^l)] ~P cos(cpr-cpp) ^^^^ ^/ 1+P»^ 

cos^'da= + PdqPpdS J i^yts 

and, with I - length of body, 

, ^^"^ , [r-r«(x-e)] -P cos(<Pr-.q)p) 

2Trf (x,r,<5Pr)+ // f^.P.cPp) ^~-3 ^ 

o o ^ 



VlW^ „ ^ , ■■ J TT cos cp 

P d <Pp d 5 = ^ ^.^ (26) 



whicli, abbreviated, yields: 

(x - 1)^ + + P= = a^1 

.+ 2 r P = i^^, 

<?'r •* 9p = X, . 



(27) 
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Bo that 

= (x - i)^ + r3, + p2 - 2 r P cosCcpj, - 9p) 
== 1- c6B(qpp - 9p) = as - gob. X.. 

It is ,r.e^dily seen that f(x,r,q5j.) must "be propor- 
tioaal to W cos cpj.. Consequently, with 

f(x,r,q5j.) = f®(x,r) W cos cp,,, ' 

a few changes result in 

P ^ ^""^ e r r . ^ cos X ■. - 

2Trf®(x.r)+ // f^d.P) |r-r'(x-|) I 



0 O 
CU S^ X, 

- P 



a^- cos 



cos A. -I 



ilEJL p dXdl _ ^ (28) 
l+r»^ Vl+r«^ 



In this integral equation the integration can be eS:« 
tended over X . It is necessary to define 

Jx = T — ^SlA cLX (29) 

° 'v^ - GOSX.3 



and 



2TT cos'' X ■ ^ , ^ 

Js = / d X . (30) 

0 ^Aa^ - cos X^ 



These integrations lead, after a numlier of transformations, 
to elliptical integrals. Putting (reference 4) 

X ^l-t^ 2dt 

tan - = t, cos A = ; — TTTf °- - ' — r~rs" 

2 1 + t 1 + t 



and ahhreviating. 



1-4 
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yields 



4 



OO 



d t 



8 a' 



U3 

/ o 



d t 



4 a* ^ 

+ ■ ■ . , ■ — 5 — / 



d t 4yas + K>2 dt/^F+t^ 



K* 0 



Moreover, the introduction of 



t^ = x: 



1 - T' 



and = 1 - \^ 



results, when K denotes the complete elliptic normal in- 
tegral of the first category, and E that of the second. 



J, = 



Js = -r Ji + 

.,8 



K*y a^+K^ 



For the sake of simplicity, we put 
4 



Ji = 



K^y a^+K^ 



(29a) 



[&i (ks) + asCk^)] . 



(3Qa) 
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that is. 



(j^(k*) = -~S!_ .E(k«) -E(k2)= i §::^ .E(k^)-K(k^.) (31) 
aa-KS 2 1-k^ 



2 k*^ ■ 
■ = 2 -pr- ■{E(k^) - E(k^)j. + ^E(k«)+K(k^)j. (32) 



Then (28) 'becomes 

I , . 

2TTf^(x,r)+ 1/ f® (e,P)-7==ii==^-fGi(k^)r-r« (x-e) + (r-P)l -PG2(k^)) 



l+r' V-i+r« " 



In this manner the two-dimensional integral equation (22) 
has been reduced to a unidimensi onal one, with nucleus or 
core: 

K(x,rje,P) =^ i&i(k^)[^r' (3:~e) + (r~P)3 -. 



This core K is examined more closelyi" It contains 
the^fun^ions Gri(k^) and G-a (k^) , which are aggregates 
of wlT&io elliptic normal integrals. Both hecome infinite 
for k^ 1 , Therefore, it is necessary to examine the 
behavior "of the core in the vicinity v'Ie .=.1- .. more close- 
ly. Conformably to the assumptions made previously, a reg- 
ular fixed point x,r is stipulated (points x = 0, r = 0 
and X = I, r = 0 are excluded for the time being). The 
term 



lT;AaC!.-i^.. ieoiinical Mempipanduin^lTo', , S75 



goes, when denoting (| - x) by |®, toward 1 as 
1 - 1®^, where - c^ ~ ( ^4 ' ^i*g ' ^° long as 

4 r^ ^iSax* !Clie term [- r« (x - |) + (r - P)] goes to- 
ward 0 as - ~ I®",,, 

Since K(k^) for k^— » 1 goes infinite (reference 



5) as - 



E(3j:°) = 1] , it is 

limit Gi(k=) [ -r' (x-e) + (r-P)] = ^j-^ = 1 
. s , 4 C ■ 1+r 

and . I. (35) 

limit P Gs (k^) = limit P 3+ln 

Thus, 'core Z(z,r;^,P) goes for k* — *1 or |® — 0 

' \JL 

. 4 ■ 

For the points x = 0, r = 0 and x = I, r = 0 k' ( r) 
"is = 0, as far as point t = 0 or \ - \ itself, where 
k^ = 1; that is, K(x,r;|,P) = 0 for all points outside 
of the fixed point. Consequently, the integral 

I 

J Z(x,r;1,P)f® d| 

0 

vanishes for tliese points. The stern of an airship is us\i~- 
ally pointed, obviating the necessity of deliberating 
whether the limitation imposed previously can be removed; 
but, as a matter of fact, it cannot, (Eeference <j) Prom 
the practical point of view this is of no significance, 
•since the stern of an air ship has .no .mathematically exact 
tip, but raer.eiy a point of very pronounced but withal fi- 
nite curvature* 

-. ;Thu.s , abbreviated, .as on the preceding page, the in- 
tegral equation becomes 

I 

2TTf®(x,r)+ / f®(^,P)K(x,r; e,P)de = — (35) 

0 ^/ 1+r '2 

*The s-abscript "max" is to indicate ;the maximum for ^ 
at which the interpolation formula is still valid. 



In 4 i» ^-remains finite, and [limit 



to infinity as In .' '. - • that is, it is integrable. 
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To resolve" it, the well-known iteratioia method is expedi- 
ent, Oho.ose as 0th approximation 



2 TTyi':+ !P:<3 



•q . I- . . . . . 

define with this, ^ f q the integral 

f ® Z d g / f ® K d I • 
0 0 

and thus compute the first approximation f®. .The ihser- 
tioij.- of this value into the integral yields the second ap- 
pros:imati on , etc. The convergence of ■ the method is as- 
sured (reference 7), because fo(x,r) is f inite .'''Nucleus 
Z(x,r;|,P) is logarithmically infinite at one point only 
and is less than- 1 ; f ® (x, r) _i^ given "by a series of 



rising powers of i/tt. I-li' Ijuttoiuu-w 



E(x,r;e,P) = — 7.4== /T+P^.K®(x,r;|,P)- ■ ^tS-S)- 
r/l+r''^ ■ . V 



and 



(x.r).= / K®(x,r;e.'p)de (39) 

0 ^ 

is initially put*, with 

Pi.(.x.r) =/ K®(x,r;|,P)di . (40) 



fe(x.r)=f«(x,r) --5 



Pi (x,r) Pg (x,r) P3 (x,r) 



TT TT^- Tf° 



(41) 



*S'or r = 0» see preceding page. 



3^8 
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; 3, . -The Practical - Solution of the Integral Equation 

a) Numerical definition of core function.^-' The most 
time-consuming part of the solution is the definition of 
the nucleus or core K^(3:,r;^,P) for different fixed- 
points (x.r) of the "body. The functions (k^ ) and 
C^sC^^) heing aggregates of whole elliptic normal inte- 
grals, are so plotted against , as to "bo accessihle. 
It should ho romemberdd that "both functions^oo for — 1 . 
G-3'(k^')' ■bocomes logarithmically infinite, that is, compar- 
ativoly inferior to Gr^ (k^) , which approaches infinity as 

— .1 — Por moro. oxact reading, it is therefore advisahlo 
1 k 

to plot tho .function H (l:^ ) = (l - k^ ) (k^ ) for Is? val- 
ues contiguous to 1, which romains finite for Ic^ — >-l, 
Thon mark off H(k^) for a prodetormined k^ and divide 
this, value hy (l' - k^,) , which yields 

Figures 10, 11, and 12 show the course of Gi (k^) , H(k^) 
and GsCk®). The calculus for a predetermined "body then 
"begins with' the det ermina-ti on of the "form factor" k^ (x, 
r;|,P) for a fixed point (x,r) (fig. 13): 

kMxr-e p)=i-5i!r2Sl=i-i"rJ2.J:J5iiPZ=i-li (42) 

k (x.r.^,P; ^^^^^ ^^^^^ ^ ^^_^p^ 1 

k^ can "be computed "by measuring off the distances 
rQ and r-j^. It was found, however, that the purely nu- 
merical determination takes about the same time and yields, 
of course, more accurate data, especially in proximity of 
(x,r). 



*k is then written in the form of 



1 + 



/r + 



N.A.C.A. Technical Memoranduni •Sro.--675 



15 



^ , Without mathematical terra given ' for ■ the contour, the 
calculation df 



QCx.r-.e.P) .= [ (r Py - r»(ac I)],. 



offers serious obstacles. 

This term Q stipulates very exact knowledge of- the 

Ay* 

tangential direction ^' = J^' (fig. 14) , and of the sec- 
ond derivation r", in order to define the limiting value 
q G-i (k^) for kf — *-l. In Figure 14 Q corresponds, to 
distance AP. It is readily apparent to 'What errors any 
inaccuracy in r' of say, 2 per cent, can lead on a long, 
slender "body. The at solute error of Q grows proportion- 
ally to (x - I), Of advantage is it that is very' 
small (Gi (k^) also) for large (x - ^) .. In close prosiai- 
ity of X, however, ~ that is, for small (x - |) , G-i 
(k®) ass-ames very high values, so that relatively minor 
discrepancies in have a .pronounced effect- on 



■because of the existing difference. If the meridian curve 
of the body is given empirically, a very careful determi- 
nation of the first and second derivative of r is abso- 
lutely necessary.. In a given case, it is necessary to 
again integrate* r" and r* 'obtained by numerical dif- 
ferentiation, to ascertain whether the initial values 
* r(x) are attained again - within the permissible, limit s. 
If it results in an inadmissible discrepancy ArCx), it 
is possible in most oases to approximate it by- a simple 
algebraic function ajid in this manner effect the necessa- 
ry corrections on ir" and r'.. 

b) Treatment of point of discontinuity of the core.- 
After computing the core for a number of fixed points 
(x,r) - the integral can be determined. The core 

being logarithmically infinite for i — x or k^ — »- 1 , 
integration is possible in interval also, which con- 

tains ■ X, provided . f remains finite. ■ Whereas the in- 
tegration can be made graphically (or planimetrically) 



K(x,r;e.P) = I Q— P G? 




♦The integration of empirical functions is admittedly much 
safer than the differentiations as calculation method. 
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outside of thi s' interval j within this critical interval 
it must te achieved numerically. (Pigure 15,) The plot 
shows the curve for K(x,r;J,P), Assume the graphic inte- 
gration is made from 0 to x - A and from x + AXj. 
to .1, The choice need not necessarily "be Axr = Ay^i , 

hut for expediency, we assume it is, and express it as 
A X, Then the integral 

x+Ax 

1= / ^^ii,P)L^(:^>Tii,9)ii (43) 
x-Ax ■ 

must he computed. Herein 

P. (I I P) =1 for computing P]_ (x,r) ■ 

■ ■ • . " ■ (44) 

P^(e,P)= ^^^' « " ife (x,r) etc.. 

lov its value is written: 

x+Ax -ft (^-.P) , . . 
1 =- / - 7 - jQi Ck^) [(r-P) -r'(x-.|)] -P Gs 1 d I .. 
x-Ax V a^ + K^ . -» 

For I — >-x, ' . ' 

Gi(L^) .[ (r ^ P) - r« (x - D] 

remains finite, as shown elsewhere thus integral 
,x+Ax 

Li 7 -TT—rr Pn(e,P)Gi(3^^) C (r-P)-r' (x-|)] d e " (45) 
x-As Vcc^+K^ 

can he graphically determined* There remains 
x+Ax 

Ls = / Pn * P) P-Cfe(3g^) a I ■ (46) 

x-Ax y a^-K^ 



Since G-s (k®) approaches infinity as In 
■ ■ o -4 
when — > 1 , we put 
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The .function 

K(k3)+ln 





4-3 


■ 4 





(.47) 



remains finite for 'Is?-^ L,.';'; (Figure 16,) 
Integral Ls resolves itself to 

lis ~ Isi + LsS +'I'33; 

W'lt.h ■ " 

x+Ax • ■ ■ 

Ls 1 = / = Pn(^.P)P 1^(1=^)* ^ C48) 

x-Ax y a2+K2 

* 

4 r P 

and, taking into account that, k® = — r 

+. 

x4Ax 

L,3 f Cp^(|.P)-p^(x.r)] y?-f ±±-^-.in.^4=^;d ! 

x+Ax , 

, . 1,33 = / ) In d e (50) 

2y r x-Ax ^ " 4 

Lsi and Las are to te defined graphically, The 
integrant of Lag remains finite, provided Pj^dtP) is 
a continuous function, For expedient evaluation function 

M(3^) = m A^"^ ' (51^ 

. ■ . ■ . . k . 4 . . 

is constructed. (Pig. 17.) 

Integral Laa is numerically interpreted. Develop- 
ing 

=1 - +< 7 e^^ + A 

,e3 se3 .1- 

P - = r...+ r« + r" 4r- ^'"-V- + J (52) 

2 v . 6. . 



♦Follows from easy extrapolation of equation (42) 



S2 
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yields 
where 



4 - 3 k' 



(53) 



ft - 3li 
Pi ~ 2r 



ps = 3,5 e' 



4r' 8 V r y 



Pa = 3,5 



2r 



-37 



(54) 



As a rule we stop after tjxe term with. and use 

the term of the 3d order to estimate the error. Stopping 
after the terms of the 2d order manifests for L33 (50) 



4 



A xpjj (x-.r) 



(1+ Ax2)ln ~ +(ln Ax-1) 



+ ^- Ax^ (m Ax - |)1 

Tile inclusion of the terms of the 3d order in (52) would 
yield a correction of the order of 



Co = Axpn(x,r) ^ (^^ 



As? , o Ax^vf 



In this manner a reference point is olstained for the se- 
lection of A'X,. 

c ) Examples ,- The development of the formulas for com- 
puting the surface, superposition ws,s followed ty various 
illustrative examples in order to check the utility of the 
method: • ■ . 

1, By a solid-of rotation produced at distance 2rma,s 
"by a source and a sink. 



*Por computing this correction, we put 
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. 3,, By. .Puhrmann "bpcLy Ko, IV* 

3, By Piihrmann "body Uo," 1, 

Tlae last was chosen tecause of its specially inconvenient 
aspect when computing the axial superposition,, 

The different approximations of the function of super- 
position f®(2,r) have been reproduced in Figures 18-20, 
The convergence of the method is seen to "be quite good. 
In example 1, the single correction has always the same 
prefix throughout the entire hody length, so that the er- 
ror caused hy the alternating prefixes of the corrective 
terms (41) is readily estimated. The first correction,. 
of course, produces the principal change in all 'bo'dios. 
Compared to the first calculation of the .core, the deter- 
mination of the higher approximations presents few ohsta- 
cles, For examples 2 and 3, the corrections change the 
prefix along the body - an effect of the almost conical 
stern. This change in prefixes is by. itself unfavorable 
for the convergence-. However, it being superior for slon~ 
der bodies to that of example 1, the significance of change 
in prefix is withal subordinate* 



0. Calculation of Flow 



a) Tangential vel o city on the surfa;ce pf the body , -■ 
In practice it is important to be able to compute the 
pressure on the surface of the body, which meians that the 
velocity must be knoWn at every point of the surface. The 
velocity is divided into its components normal and tan-^ 
gential, .The first, induced by the superposition of the 
parallel flow in direction of negative axis z over the 
source flow, is zero* The second .i s divided in a point 
(x,r,q)) into components % and Wg at right angles to 
one anothero Component w^^ clings to the parallel of lat- 
itude in direction of the tangent, component Wg in direc- 
tion of the tangent of the meridian through point (x,r,<p) * 
(Fig, 21,) ■ is figured positive in direction of grow- 

ing angle qj ; Wg • is ..valid as . positive, when pointing from 
nose to stern* Now, let f®(^,P) denote, as before, the 
surface superposition conjugated to velocity W = 1, It 
results for w^ and Wg in' 
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• - • I F , » . 

TTx « W Sin qj [1 + — — ^h=rr / f^CI.P)*/! + P'^- 

X --— i:^:^ &3(k^) dl] . ■ (55) 

ya^ + 

m^lxere 

■ %(^^) = {(2 - k^) KCk^) - 2 ECk^Xj. (56) 

(fig. 22) and 

* ■ ■ r "p ' ? 

Wg = W coscp J i — + — 

■ "■ V. V 1 + r«3 ry 1 + r«2 

o >y + k2 ■ . ■ > 

X Gi(k^) - P r« Gs(k2)] dl| (57) 

G-3(k^) is 1 o gar itlimi.c ally infinite, the integral in 

(55) therefore finite, Hote, when evaluating the integral 

o'f (57) that, whereas Gg (k^) . is logarithmically infinite, 

Gi(k^) is as* — ^ = — — — r . On the other hand, since 

. . 1 - k= €^ ■ 

<f = (x-l) + rt(r-P) = ~ (l+rt^>p ^ ^ + ... 

changes its prefix for |® = 0, Q®' tJ^ (k^) "goes, when 
proceeding from the left on x,. toward + °° , and toward 
- CB when approaching x from the right, Ihe integral 



* ■ 

,mua.t therefore be expressed hy its "principal value" 



* f>e - X. 
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limit y tsLsHLxs^ ^ di 

Another point to "be noted with Wg is that the inte- 
..gral for points x = 0, r = 0- and x = I', r = 0 does 
not disappear (in contrast to the normal component set up 
hy source superposition alone). It is hest to define the 
integral for this point . separately , foregoing the reduc- 
tions, which lose their significance because is al- 
ways zero up to one point. Tangential v^elocity % set 
up in point., x = 0, r = 0 hy source superposition "is in 
direction of negative axis z and known in magnitude hy 
(fig. 23) 

_ , +Tr/2 I f^(e,P)coscpp ■ > , 

Wgl = 2W / / g =- sin ^yi+P'^dl cos(pp Pdcp, . 



-TT/2 o 

= 2¥ 



-Kr/2 



} fe(i,p)p^y i+p»^ , t ; 2 . 

y 3^ d i S cos^ (pp dqp 

o ^ -17/2 



= TTW / £!Ck£l£!Zi±Pl!M . (58) 



wherein H = + P^ . 

b) The potential in any point of the space .- This is 
given hy 

$ = - 4 W cos (ft / ^^<^' P> '^'^■^ 211 G3(k^)P dl (59) 

0 J~(l^ + K2 . . 

G-3(k^) heing the function defined by (56). 

/ r \^ 

JTote, that, = 1 - ( --2-) is always ' 1, T7hen com- 

puting the potential in points outside of the body. As e, 
result, no discontinuity places of the integrand appear. 



2S 
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An outline of Von, Zarman' s metliod of computing the 
potential flow' of airships in yaw "by means of partially 
constant dipolar superposition on .the axis of the "body is 
followed by several considerations for.. "beginning and end 
of the superposition. Then- this method i s improved hy 
postulating a continuous, in part linearly variable di~ 
polar superposition on the axis. 

The second' main part of the' report brings the calcu- 
lation of the. potential f low by mean.s of sources and sinlcs, 
arranged oh the surface of the airship body. The integral 
equation which ^gu^ satisfy this surface superposition is 
posed, and • the .fojao, reduced to functions developed from 
whole elliptic normal integrals. The functions are shown 
diagrammatically . The integration is resolvable by iter- 
ation. The consequence of the method is good. The for- 
mulas for computing the velocity on the surface and of the 
potential for stnypolnt conclude tho report. 

Translation by J. jranier.,,,. 
ITational Advisory Committee 
for Aeronautics. 
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Figs. 1.2,3,4 




P(x,r) 




Fig, 1 Coordinate axes 



Fig, 3 notations 




Fig. 2 Body with dipolar super- 
position on the axes. 




Fig. 4 Intervals "by partially 

linearly variable di- 
polar superposition. 
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Figs. Say^fea^eiStec 




Pig. 5a 




Fig. 5D 




Fig. 6a 




Figs. 5,6 

Computed 

dipola,r 

superpositions 



Fig. 6c 
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Pigs. 7,8,8a,8TD,8c 
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Pig. 8 
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Pig. 8a Pig. 8l3 Pig. 8c. 

Pigs. 7,8 Compated dipolar superpositions. 
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?igs. 9,10,11 
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Pigs. 12,13,14 




■.rig. 12 The GpCk^) fumction. 




Fig. 13 Definition of Pig. 14 n-eomeiirical signif- 

forin factor icance of Q 
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Figs. 15,15,17,18 
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Fig. It) Jxmction of curve K(x,r,|iP) '^^^ 

for x= 50. (Fuhrmaim shape 
Mo.l, 1= 230) 
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ig. 16 The N(k2) function. 
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Eig, 18 Siirface superposition of 
source - sink "bocly. 
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Pigs. 19,20 
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Pig, 19 Surface soperpositioii of Pulinaoiua shape Fo. IV 
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Pig, 20 Surface superposition of Puhrmaim shape ITo. I 
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Tigs. 21,22,23 




Pig. 21 Velocity components on surface of "body. 




Fig. 22 The Gr,(k2) ftmction 




Pig. 23 Meridian cut at angle cp 
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